Scattering of Graphene plasmons by defects in the 

graphene sheet. 

Juan L. Garcia-Pomar,* i * Alexey Yu. Nikitin,^* and Luis Martin-Moreno*'*'* 

Instituto de Ciencia de Materiales de Aragon and Departamento de Fisica de la Materia 
Condensada, CSIC-Universidad de Zaragoza, E-50009, Zaragoza, Spain 

E-mail: lmm@unizares] 
Abstract 

A theoretical study is presented on the scattering of graphene surface plasmons by de- 
fects in the graphene sheet they propagate in. These defects can be either natural (as domain 
boundaries, ripples and cracks, among others) or induced by an external gate. The scattering 
is shown to be governed by an integral equation, derived from a plane wave expansion of the 
fields, which in general must be solved numerically but it provides useful analytical results for 
small defects. Two main cases are considered: smooth variations of the graphene conductivity 
(characterized by a Gaussian conductivity profile) and sharp variations (represented by islands 
with different conductivity). In general, reflection largely dominates over radiation out of the 
graphene sheet. However, in the case of sharply defined conductivity islands there are some 
values of island size and frequency where the reflectance vanishes and, correspondingly, the 
radiation out of plane is the main scattering process. For smooth defects, the reflectance spec- 
tra present a single maximum at the condition k p a « y/2, where k p is the GSP wavevector and 
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a the spatial width of the defect. In contrast, the reflectance spectra of sharp defects present 
periodic oscillations with period k f p a, where k' p is the GSP wavelength inside the defect. Fi- 
nally, the case of cracks (gaps in the graphene conductivity) is considered, showing that the 
reflectance is practically unity for gap widths larger than one tenth of the GSP wavelength. 

Introduction 

In the last few years, it has become evident that graphene not only displays remarkable electronic 
properties but can also play a significant role in photonics.^ One aspect that has recently attracted 
much interest is that doped graphene supports bound electromagnetic modes, known as surface 
plasmons (GSPs), 3 that have the appealing characteristics of being both confined in a length scale 
much smaller than the free space wavelength,EH3 and potentially controllable using external gates. 
Very recently, the existence of highly confined GSPs has received experimental confirmation.^^ 
Several aspects of GSP have already been studied theoretically, such as the efficient and direc- 
tional coupling with nano emitters (and the associated enhanced spontaneous emission rate),^^ 
enhanced absorption and resonance diffraction,^^]] metamaterials and antenna applications'^!!! 
and their wave guiding capabilities in 

ribbons'2M 

and edges. 1 1 

However, very little is known about how GSPs behave when they encounter defects in the 
graphene sheet they propagate in. These defects can occur both (i) naturally as, for instance, 
kinks appearing due to fabrication process,^ domain borders in graphene growth by CVD^the 
presence of multilayer islands, 32 cracks^^Sl anc [ different domains in CVD graphene^ or (ii) 
be created externally, for example as the changes of conductivity in gate-induced p-n^ or p-n-p 
junctions.! 37 * 3 ^ 

In this article, we present a theoretical study of GSP scattering by one-dimensional conductiv- 
ity inhomogeneities. Calculations are conducted with an original method based on the Rayleigh 
expansion, which has the advantage of providing analytical expressions in some limiting cases. 
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Figure 1: (Color Online) (a) Schematic geometry of free-standing graphene with a defect zone 
(blue color) that interacts with a GSR The figure shows the conductivity profile (bottom) and a 
representative case of the computed the total (incident + scattered) magnetic field modulus \H\ 
(top), (b) Diagram for the different types of conductivity defects analyzed, (c) Schematics of the 
scattering processes occurring when a conductivity defect is present in a graphene sheet. 

We consider a free-standing graphene monolayer, placed at z = 0, with a spatial inhomogeneity 
in the two-dimensional conductivity a. Actually, the presence of a substrate may be indispensable 
for applications, but it does not change any of the fundamental scattering properties of GSPs (af- 
fecting mainly the mobility of the charge carriers), which is why in this paper we concentrate on 
the simplest structure. We will analyze one-dimensional (ID) inhomogeneities, with translational 
symmetry in the in-plane direction perpendicular to GSP incidence (the j-direction), so a — a{x). 
The geometry of the system is schematically shown in Fig. [TJ 

Away from the defect the (frequency-dependent) conductivity of graphene is Oq. Defects have 
a characteristic width a, and a conductivity o m at the defect centre (x = 0) which can, alternatively, 
be described by the relative change in conductivity 8 = (<7 m — Gg) / Og- 

In the scattering geometry, a monochromatic GSP (time dependency e~ lC0t ), propagating along 
the Ox axis (from the region x < 0), impinges the defect, which induces some reflection back into 
the SPP channel, as well as some radiation out of the graphene sheet. Notice that, due to the 
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symmetry of the problem, all scattered waves have the same polarization as the GSP (transverse 
magnetic). The scattering amplitudes can be computed by using numerical solvers of Maxwell 
equations. Here, we present an alternative method, based on the Rayleigh plane-wave expansion 
(RP WE) ,122140] w ] 1 i c ] 1 i n general must also be solved numerically but presents the advantage of 
providing analytical expressions for the scattering coefficients in some limiting cases. We leave 
all derivations for the Supporting Info and present here the main equations. Within the RPWE 
method, the electromagnetic field is written (all other components can be readily obtained from 
Maxwell equations) as: 

/oo 
G(q)B(q)e^ x dq (1) 
-oo 

where q and q p are x-components of the wavevectors of a plane wave and GSP respectively, nor- 
malized to the wavevector in vacuum g — 2%jX. B(q) are scattering amplitudes, which satisfy the 
integral equation: 

/oo 
Aa(q-q')G(q')B(q')dq' (2) 
-oo 

Here Aa(q) is the scattering potential related to the Fourier transform of the inhomogeneity of the 
dimensionless conductivity Aoc(x) = (2k/c)(g(x) — Gq) and G(q) = q z j (l+q z ac), where c is the 
speed of light and q z = y/l — q 2 . 

The GSP reflectance and transmittance (R and T , respectively), as well as the fraction of energy 
flux scattered out of plane S can be obtained from the amplitudes B(q) (see Supporting info). For 
instance, R = \2nB(—q p ) / (qpOCo) | 2 . 

It must be pointed out that this model does, in principle, account for losses in the graphene 
sheet. However, for the frequencies and defects that will be analyzed in this paper, for which a is 
much smaller than the GSP absorption length L p , the inclusion of losses leaves the scattering coef- 
ficients virtually unaltered. In a real situation in order to obtain reflexion/transmission coefficients, 
one has to normalize the GSP amplitudes by the decay factor e~ d / L p, where d is the distance that 
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GSP runs from the launching point to the destination where the amplitude of the scattered GSP is 
measured. We have checked (running simulations with both the RPWE method and a commercial 
finite-elements code 41 ) that once this procedure is applied for the lossy case, the scattering coeffi- 
cients virtually coincide with those obtained when losses are neglected. Therefore, and in order to 
concentrate on the scattering coefficients intrinsically due to the defect, all calculations presented 
in this paper have been obtained setting Re[<7G] = 0. In this way, current conservation implies 
7? + T + S = 1. Throughout the paper the conductivity is taken from the RPA expression^HSI anc [ ? 
for definiteness, we consider that the chemical potential is /I = 0.2eV. As we will show, this choice 
is not essential, as most results only depend on ji through the value of the GSP wavevector. 

Of course, for the scattering coefficients associated to any particular defect will depend on its 
conductivity profile. Here we do not attempt to computing this profile; instead, we will assume 
some basic spatial dependences for the conductivity and compute how they scatter GSPs. We 
analyze two differentiated main cases: (i) smooth variations in the conductivity, described by a 
Gaussian profile a{x) = <7g{1 + 8 exp(— 4x 2 /a 2 )}, where a is the full spatial width at l/e relative 
conductivity change, and (ii) abrupt ones, represented by a step defect a(x) = cjg{1 + 8 &(a/2 — 
|*|)}, where 0(jc) is the Heaviside step function. 

Smooth defects 

Let us advance that for, all smooth defects considered, our calculations show that the scattering 
out-of-plane is extremely small (S/R< 10 -4 — 10 -2 ). This point will be discussed later on; now it 
allows us to focus on the reflectance, from where the transmittance can be obtained as T « 1 — R. 

We first consider a Gaussian profile variation in the conductivity with a small <5. This profile 
provides a good approximation to the realistic variations on graphene conductivity arising from 
atomic steps in the substrate, when graphene is grown on SiC.^ In this case the characteristic 
width of the Gaussian profile is a ~ 20 — 50nm, and the conductivity relative change 8 ranges 
between —0.01 and —0.5. 
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Figure 2: Reflectance spectra R(X) for a GSP impinging a shallow Gaussian conductivity defect, 
for different defect widths and two values of <5. Inset: R as a function of a/X p for 8 — —0.2 and 
different widths, showing that the reflectance scales with a/X p . The results obtained with the finite 
element method (open circles) for a — 30nm and 8 = —0.2 are also shown, in order to confirm the 
validity of the RPWE calculations. 

Fig. [2] renders the reflectance spectra R for defects with different widths and two values of the 
relative change in conductivity. We observe that R has a maximum for each value of 8 and a. This 
maximum arises as a compromise between two distinct asymptotic dependencies. For small GSP 
wavelengths, X p , the GSP follows adiabatically the variation of the conductivity and virtually no 
reflection is generated. Conversely, in the long-wavelength region, the defect width is very small 
in relative terms (a <C X p ) and so is the reflectance, which decreases with X due to the decrease in 
a/X p . In between these asymptotic decays there is maximum for 7?, related to the g-space Fourier 
image of the Gaussian conductivity profile. 

In order to gain more insight into the behavior of the GSP reflectance, and obtain some quan- 
titative estimations to support the exact calculation, we compute the plane- wave amplitudes B(q) 
within the first-order Born approximation (FOB A). This approximation is valid for small variations 
of 8 and corresponds to neglecting the integral term in the right-hand side of Eq.(2), keeping only 
the linear term in Aa. 

Within FOB A the scattering amplitude reads B F0BA (q) = - Aa (q - q p ) . Using (Xq 2 = l-q 2 p ^ 
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Figure 3: R/8 2 as function ofa/X p , for for a Gaussian conductivity defect with a = 50nm and dif- 
ferent values of <5. We observe that, as 8 decreases, the curves converge to the universal reflectance 
spectra predicted by the FOBA. 

— q 2 we obtain (see details in Appendix) 

R FOBA = ^ {kpa) 2 e -^k p af 5 2 Q) 

Notice that, within the FOBA, the quantity R/8 2 is a universal function of a/X p . It also predicts 
that the maximum in reflectance occurs when a/X p — 1/ (V27t) « 0.22 (independent of <5), with 
a maximum reflectance R^^ A — {%/2e)8 2 « 0.58<5 2 . The validity of the scaling of the reflctance 
with a/X p predicted by the FOBA is shown in the inset to Fig. |2j for different defect widths and 
5 = —0.2. Additionally, Fig. [3j renders the scaling with 8 computed for a narrow defect, together 
with the prediction by the FOBA. The FOBA captures very accurately the spectral position of the 
maximum reflectance, even for moderate variations in conductivity, and gives a good approxima- 
tion to the full reflectance spectra. The FOBA also provides insight into the relative strength of 
reflectance and radiation channels. The scattering strength depends both on (i) the density of final 
states (which is much larger for GSPs than for radiation channels) and (ii) a matrix element, given 
by the Fourier component of the conductivity variation evaluated at the wavevector difference A& 
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between the GSP one and that of the final state (i.e., Ak = —2k p for reflectance and Ak « — k p 
for radiation processes). In the case of smooth defects, the FOBA shows that the density-of-states 
factor dominates over the "matrix element" one (see Appendix). Actually, at the reflectance maxi- 
mum, the FOBA predicts S/R = 0.5e 3 / 4 \ (Xg\ 2 , which is in the range ~ KT 4 -l(T 2 for the values 
of (Xg relevant for GSP propagation (ccg ~ 1 — 5ao, where (Xo w 1/137 is the fine structure con- 
stant). This preponderance of 7? over S holds even for larger defects strengths, where the FOBA is 
no longer strictly applicable. 

It is also interesting to analyze the scattering by smooth defects with large <5, as they can be 
produced by changing the carrier concentration in graphene (and thus the conductivity) with an 
external gate J 37 * 38 ! The Gaussian shape in this case may simulate an~ —n — n~ (or p + — p — p + ) 
junction. In this case, the defect width will depend on the geometrical details of the gate, but can 
be expected to be of order of 0.2 — 1 jim, and 8 can be considered a tunable parameter ranging 
from — 1 to 0. 

Fig. [4] presents the reflectance spectra for different values of <5, for the fixed defect width 
a = 400nm. These results show that, for large relative changes of the conductivity, the maximum 
reflectance occurs for even smaller defect widths than those predicted by the FOBA and that, for 
large 1 8 | , there are spectral regions where the reflectivity is high. This point is even more apparent 
in Fig. [5][a), which renders the reflectance spectra for 8 = —0.9 and different defect widths. Inter- 
estingly, the reflectance still satisfies approximately the scaling relation in ajX p predicted by the 
FOBA, see Fig. [5^b), although for this large value of |5| the FOBA is no longer a good approxi- 
mation to the scattering amplitudes, which must be obtained by solving the full integral equation 
Eq.(2). 

It is remarkable that for the large relative change in conductivity considered in Fig. [5} the 
reflectance in the spectral region ajX p > 0.2 is small. Actually, there are values ofa/X p where the 
reflectance vanishes and, given the scattering out of plane is negligible, the transmittance is almost 
unity. As the GSP extension in the direction perpendicular to the graphene sheet scales with the 
conductivity, the GSP is very strongly bound at the defect centre. Then, unit transmittance and 
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Figure 4: Reflectance spectra for a GSP impinging a Gaussian conductivity defect for different 
values of <5. The defect width is a = 400nm. 

conservation of energy imply that the electric field is strongly enhanced at the centre, as illustrated 
in Fig. |6]^a). The scaling of the electric field amplitude can be obtained by assuming that, for 
adiabatic propagation of GSPs, the electrical current along the graphene sheet is constant, i.e. 
\J(x) | = | a{x) | • \E x {x) | = constant, leading to \E x {x) | 1/| o{x) \ . As for the other EM components 
of the GSP field, we know that they satisfy \E X \ ~ \E Z \ ~ \q p H\. So, taking into account that locally 
q p (x) ~ i/a(x), we arrive at \H(x) \ ~ \a(x)E x (x) \ oc constant. The spatial dependence of the GSP 
field is rendered in Fig. [6fb), together with the conductivity profile, fully confirming the predicted 
scaling behavior. 



Abrupt defects 

One paradigmatic case of defect with abrupt change in conductivity is an island of multilayer 
graphene, placed on a graphene monolayer. As for small number of layers the thickness of the 
multilayer is much smaller than the GSP extension along the normal to the sheet, this thickness 
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Figure 5: (a) Reflectance spectra for a GSP impinging a Gaussian conductivity defect, for different 
widths and for both 8 = —0.9 and 8 = —0.4 . (b) R/8 2 in function of ajX p for the same pa- 
rameters considered in (a), showing the approximate scaling of the reflectance spectra. The open 
circles render calculations performed within the FEM (in order to validate those performed with 
the RPWE method), for a = 400wn and 8 = -0.9. 

can be neglected and the multilayer region can be approximated by a conductivity defect. Here we 
analyze the scattering by both bi- and tri- layer strips, and approximate their conductivities as twice 
(8 = 1) or three times (8 = 2) the conductivity of a monolayer, respectively. The FOBA calcula- 
tion for these rectangular-type defects gives R F0BA — ^iv?{k p a) <5 2 , predicting that the reflectance 
spectra oscillates periodically when expressed as a function ofa/X p . Within the FOBA, reflectance 
minima occur at a = nX p /2, n— 1,2,... (which can be interpreted as the constructive interference 
in the backward direction between GSP partially reflected at the edges of the defect). The FOBA 
is not a good approximation for these defects where 8 is not small, as it fails to take into account 
the modification in the GSP field inside the island due to the change in conductivity. However, we 
have found that the full calculations follow quite approximately the periodic behavior predicted by 
the FOBA, but as a function of a/X' p9 where 7J is the wavelength of the plasmon corresponding 
to the conductivity inside the defect, see Fig. [Tfa). Notice also that, for these abrupt defects, the 
scattering out of plane also presents an oscillatory behavior. Its amplitude, although smaller than 
the one for the reflectance, is not negligible, and is peaked (as also does the transmittance) at the 
spectral positions where the reflectance is minimum. 

As illustrated in Fig. |7Jb), the strong reflection of GSP at the island boundaries result is the 
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Figure 6: (a) Snapshot of the electric field norm along the GSP interacting with a Gaussian con- 
ductivity defect characterized by 8 = —0.9, a = 200nm and A = 200/im. (b) Crosscut at z=0 of 
the previous snapshot. Open squares represent the normalized magnetic field norm H/Hq (open 
squares), which is approximately unity along the GSP propagation, open blue circles render the 
normalized electric field norm E/Eq, showing the exact scaling with the inverse of normalized 
dimensionless conductivity (cc/ccg)~ 1 (continuous black line). 

formation of standing waves in the island, with a number of nodes determined by both the island 
size and the GSP wavelength there. Notice that, as the conductivity inside the multilayer island is 
larger than in the monolayer, the GSP is less strongly bound to the graphene sheet. Also, Fig. [Tfb) 
clearly shows that the out-of-plane radiation is generated at the island boundaries. 

Finally, as another paradigmatic case of abrupt defect, we study the scattering of a GSP by a 
crack in the graphene layer. The exact spatial dependence of the conductivity near the graphene 
edge is a question still under debate, with the microscopic details of the graphene edge (whether has 
a zigzag or armchair configuration) possibly playing an important roleJ 45 * 46 * Here we will simply 
assume that the conductivity vanishes within the gap region, in order to provide an estimation of the 
distances that GSP can tunnel through. Fig. [8]renders the dependence of the computed reflectance 
with gap width, for several frequencies. The results are presented as function of a/X p , showing 
that the reflectance approximately follows a scaling behavior. Still, the reflectivity is high already 
for small values of a/k p , demonstrating the extreme sensitivity of GSPs to the presence of cracks 
in the graphene sheet. The insets to Fig. [8] show snapshots of the magnetic field, illustrating (i) the 
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Figure 7: (a) Reflectance 7?, transmittance T and fraction of energy radiated out of plane S versus 
a/Xp for a bilayer (left) and trilayer (right) for A = \0jim. (b) Snapshots of the magnetic field 
norm for a GSP propagating in a graphene monolayer and impinging onto either a bilayer (left) or 
a trilayer (right) strip. The strip width is a — 1/im. 
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standing wave arising from the reflection of the GSP at the crack, (ii) the smallness of the radiation 
out of the graphene sheet (in all calculations, the fraction of GSP energy radiated out of plane is 
S « 10" 3 ), and (iii) that the GSP is fully reflected for widths a > 0.lX p . 
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Figure 8: Reflectance R for a GSP impinging onto a crack, as a function of crack width a, for 
different frequencies. FEM calculation for 6THz is shown on order to demonstrate the validity of 
the RPWE method. Inset: snapshots of the magnetic field norm, normalized to that of the incident 
magnetic field at the graphene sheet, for the cases ajX p — 0.02 and ajX p — 0.11. 

Conclusion 

We have analyzed the scattering properties of GSP by defects in the local conductivity of the 
graphene sheet. In the case of smooth spatial variations of the defect conductivity (which occur for 
instance when the defect is created by modification of the carrier concentration via a top gate) we 
have found that, for a given relative change in the conductivity at the defect centre, the reflectance 
follows approximately a universal scaling in terms of a/X p . In all cases, the reflectance reaches 
its spectral maximum value when a w Q.2X p . When a/X p is larger than that given from the pre- 
vious condition the GSP propagation can be considered as adiabatic, and thus the GSP is mainly 
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transmitted. When, additionally, the conductivity at the centre of the defect is small this leads to 
a strong electric field enhancement at the defect centre. We have also found that for these smooth 
defects, the scattering out of plane is always much smaller than the reflectance. 

A different behavior is found for the scattering of GSP by multilayer islands, placed in a mono- 
layer background. In this case, the scattering out of plane it is not negligible, although on average 
it is smaller than the reflectance. In fact, the reflectance spectra oscillates periodically as a func- 
tion of quotient between the island width and the plasmon wavelength inside the defect, with both 
transmittance and scattering out of plane presenting maxima at the reflectance minima. 

Finally we have found that conductivity gaps in the graphene sheet prevent very efficiently the 
GSP propagation, with the GSP being fully reflected for gap widths larger than ~ O.lA^. 
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Appendix 

- Derivation of the integral equation governing scattering of GSP s by a defect. 

Let us assume that the normalized conductivity of graphene a — 2kg /c is the following 
function of the coordinate x: a(x) — (Xg + Aa(x), where ccg is the background conductivity 
[ccg = oc(x —> ±00)] and Aa is a localized perturbation [Aa(x —> ±00) = 0]. The perturbation 
can be represented by its Fourier expansion 

Aa(x) = J dkAa(k)e ikx , Aa{k) = J dxAa(x)e~ ikx . (Al) 

We consider a graphene surface plasmon (GSP) propagating along a free-standing graphene 
sheet located at z = 0, from the region x < (we assume lossless graphene for the reasons men- 
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tioned in the manuscript). The GSP impinges the inhomogeneity Aa(x) localized around x — 
and is partially reflected, transmitted and scattered into a continuum of propagating and evanescent 
modes \k>. Then the longitudinal components of the electric fields above and below a monolayer 
can be exactly represented by the Rayleigh integral expansion [time dependency ~ exp(— icot) is 
supposed everywhere] 



/oo 
dkE+(k)e ik *e ik * z , z>0, 
-oo 

/oo 
dkE-(k)e ik *e- ik * z , z<0. 
-oo 



(A2) 



with k p = gyj 1 — 1/oCq, g — (o/c — 27i/X. The z-component of the wavevectors are given by k z = 
\/g 2 — k 2 (with Im(fc z ) > 0), and k zp = —gjaQ. The components of the electric field E = (E x , 0, E z ) 
and magnetic field H = (0,H y ,0) are connected through the Maxwell's equations. 

The boundary conditions at z = are provided by (i) the continuity of the parallel component 
of the electric field and (ii) the jump of the parallel component of the magnetic field across the 
graphene layer due to the induced electric current 

E+(x,0)=E-(x,0), 

(A3) 

H+ (x, 0) - H~ (x, 0) = 2a(x)E+ (x, 0) . 

The fields given by Eq.(A2) can be substituted into the boundary conditions Eq.(A3). Using 
Eq.(Al) and projecting the equations onto the basis set e lkx we arrive at the integral equation for 
the Fourier amplitude of the field E Jr {k): 

B(q) = -Aa(q -q p )- J Aa(q - q , )B(q , )G(q , )dq f (A4) 

Here we have introduced the normalized wavevector components q — kjg, q z — k z /s^ so that 
Aa(q) = gAa(k) 9 and normalized the Fourier image of the field to the Green's function G(q) = 



15 



q z /(l + a G q z ): 



E + (k) = G(q)B(q). 



(A5) 



Eq.(A4) has been solved discretizing q and replacing the infinite region of integration by in- 
creasingly larger finite limits, until convergency is achieved. A non-uniform discretization scheme 
in q-space is considered, in order to take into account the strong variations of the Green's function 
G(q)™ 

Once the Fourier image of the field is computed, we can obtain the GSP reflection, transmis- 
sion, and out-of-plane scattering efficiency. For this, we use the normalization Eq.(A5). Assuming 
that the function B(q) does not have poles in the complex plane q, the asymptotic (long-distance) 
behavior of the field at x —> ±°o is given by the contribution from the poles of G(q): 



^(x^oo, ) = (i + Tyv, 

E x (x -> -00, 0) = e ik P x + pe~ ik P x , 



(A6) 



where 



2%i 
OC^qp 



B(q P ), P 



2%i . 

-5— B(-q p ). 
CiQqp 



(A7) 



From Eq.(A6) and Eq.(A7), the reflectance and transmittance are immediately obtained: 



R 



H + tI 



1 + 



2%i , 
-^B{-q p ) 

2ni 



„3 B M 
0CQq p 



(A8) 



The computation of the energy radiated in the form of the propagating waves is done by performing 
the saddle-point integration of the fields in the far-zone, and integrating the normal component of 
the Poynting vector P = ^Re[E x H*] over a closed cylindrical surface. Then, normalizing it onto 
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the flux of the incident SPP, J p = q p \ac\ 3 /2g, we obtain 

9p\^G\ %I<1 <lz 

which can also be written in a form that permits to define the angular dependence of the radiated 
flux: 



f 7T/2 

S = / dOD(O), (A10) 

J-Tl/2 



where D(6) is the differential cross section: 

4% o 

D(0) = — : ^|5(cos0)G(cos0)| 2 , (All) 

qp\ a G\ 

with the angle 6 counted from the axis Oz. 

As we have assumed that our system is lossless, the scattering coefficients constitute the current 
conservation 

1-R-T-S = 0. (A12) 

- Scattering coefficients in the first-order Born approximation. 

The first-order Born approximation (FOB A) corresponds to the first term in the series expan- 
sion of B(q) in the perturbation amplitude Aa. FOB A can be simply obtained from Eq.(A4) by 
neglecting the integral term: 

B(q) F0BA = -Aa(q-q p ). (A13) 

The FOB A is strictly valid only for small perturbations. However, it is fully analytical and, in many 
cases, it provides very useful information. The calculation within FOBA for the defects considered 
in this manuscript are: 

- Gaussian Defects. 

In this case, a defect with width a (defined as the full spatial width at 1 je relative conductivity 
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change) has a conductivity profile Aa(x) = <5aGexp(— Ax 2 /a 2 ), where 5 = (a(x — 0) — olg)/olg. 
Then, Aa(q) = aG<3exp(— q 2 a 2 /\6) <5, with <3 = ga. According to Eq.(A8) and Eq.(A13) the 
reflectance reads 

R FOBA = | {kpa) 2 e - l 2(k P a) 2 $2 (A14) 

where we have used that (Xq 2 — \—q 2 p ^ —q 2 . The fraction of energy radiated out of plane can 
be estimated by approximating B(q) « B(0) and G(q) ~ g z for |g| < 1. With this, 

S F0BA « £ |a G | 2 (V) 2 * 4(V)2 S 2 - (A15) 
8 

Notice that the reflectance maximum occurs at k p a — y/2. Then, at the reflectance maxima the 
FOBA predicts that the ratio R/S = 2e~ 3 / 4 \a G \~ 2 « 1.8 x 10 4 /(a G /a ) 2 , where a o ~ 1/137 is 
the fine structure constant. 
- Abrupt defects. 

For a defect with constant conductivity and spatial width a, Aa{x) = 5 (Xg &(a/2 — \x\). Then 
Aa(q) = aG<3 sine (^(3/2) 5, where sinc(x) = sin(x)/x. Using this expression, the scattering 
coefficients within the FOBA can be straightforwardly computed. 
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